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Abstract 



The Lagrangian derivation of the Equations of Motion for topological static spherically 
symmetric metrics in J-(R, G)-modified gravity is presented and the related solutions are 
discussed. In particular, a new topological solution for the model J-(R, G) = R + VG is found. 
The black hole solutions and the First Law of thermodynamic are analyzed. Furthermore, 
the coupling with electromagnetic field is also considered and a Maxwell solution is derived. 

5-H . 

1 Introduction 

Recent observational data imply -against any previous belief- that the current expansion of the 
universe is accelerating Qj 2, 3]. Since this discovery, the so called Dark Energy issue has become 
the "Mystery of the Millennium" [3] and today dark energy is probably one of the most ambitious 
and tantalizing field of research because of its implications in fundamental physics. There exist 
several descriptions of the acceleration of the universe. Among them, the simplest one is the intro- 
duction of small positive Cosmological Constant in the framework of General Relativity (GR) , the 
so called ACDM model. A generalization of this simple modification of GR consists in considering 
modified gravitational theories, where some combination of curvature invariants (the Riemann 
tensor, the Weyl tensor, the Ricci tensor and so on) replaces, or is added, to the classical Hilbert- 
Einstein action of GR. The simplest class of modified theories is i^(i?)-modified gravity, in which 
the action is described by a function F(R) of the Ricci scalar R (for a review see Refs. 0E]), but 
other modifications are in principle allowed. For example, an interesting class of modified gravity 
models which may easily produce the acceleration epoch is string-inspired modified Gauss Bonnet 
gravity, the so-called i 7 '(G)-gravity jTJ [5] , where F(G) is a function of the four dimensional Gauss 
Bonnet invariant G. 

If some modified theory lies behind our universe, it is of crucial interest to extend the proprieties 
and the laws of General Relativity to its framework. In this paper, we focus our interest on the 
topological static spherically symmetric (SSS) solutions in modified gravity F(R, G)-models, that 
is, the modification of gravity is given by a function of the both, the Ricci scalar and the Gauss- 
Bonnet invariant. Typically the modified gravity models admit the de Sitter (dS) space as a 
solution, but the number of exact non-trivial (i.e., different to Schwarzshild-de Sitter one) SSS 
solutions so far known in modified gravity is extremely small, due to the complexity of the field 
equations. 
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The SSS solutions in F(i?)-gravity have been investigated in several papers (see Refs. 
EH EQ H2 US])- For the specific model F(R) = R 1+s , S being a fixed parameter, a class of 
exact SSS solutions has been presented in [T3]. For SSS solutions in Weyl gravity see also Refs. 
[T5lll6| . In Ref. [T7] a suitable derivation of the F(_R)-equations of motion based on the lagrangian 
multilpleyers has been presented and it has been shown how this approach permits to find a large 
class of solutions. Here, we would like to extend this method to the topological case of SSS solutions 
in T(R, G)-modified gravity. In Section 2 we present our Lagrangian derivation of the Equations 
of Motion (EOM) and in Sees. 3-4 we will see how this approach permits to find explicit solutions 
via reconstruction method. In particular, an exact solution with non constant Gauss-Bonnet is 
found and a large family of solutions with constant Gauss-Bonnet is discussed. 

The SSS solutions may describe the black holes (BH) and in Section 5 we will evaluate the BH 
entropy in T(R, G)-gravity via Wald method [TS] . Furthermore, it seems very difficult to address 
the analogue of the Misner-Sharp BH mass of GR to modified gravity, where the definition of a 
sensible mass parameter is a debated question. In Section 6, by following the proposal of Ref. [19| 
for .F(i?)-gravity, we will try to identified the mass with a quantity proportional to the constant 
of integration, which appears in the explicit solutions, making use of the First Law of black hole 
thermodynamics and evaluating independently the entropy via Wald method and the Hawking 
temperature via quantum mechanical methods in curved space-times [20J. An attempt to derive 
the First Law from the Equations of Motion of T(R, G)-gravity is done. 

In Section 7 we provide the formalism for topological SSS solutions in the presence of Maxwell 
field and an exact solution for a specific Gauss-Bonnet modified gravity model is found. Conclu- 
sions are given in Section 8. 

We use units offcs = c = h = 1 and denote the gravitational constant k 2 = 8itGn = 87r/Mj> ; 
with the Planck mass of M PL = G„ 1/2 = 1.2 x 10 19 GeV. 



Lagrangian approach for topological static spherically sym- 
metric vacuum solutions 



In this Section we present a suitable derivation of the J-(R, G)-equations of motion based on 
lagrangian multipleyers, which permits to deal with an ordinary differential equation system. The 
action of modified J-(R, G)-theories reads (in vacuum) 

/ <?xy/=$F(R,G), (1) 



2k 2 



M 



where g is the determinant of metric tensor, g^, Ai is the space-time manifold and T{R, G) is a 
generic function of the Ricci scalar R and the Gauss Bonnet four-dimensional invariant G, 



G — R 2 — AR^R^" + Rpy^tjR^^ 17 



(2) 



The Gauss Bonnet is a combination of the Riemann Tensor R^a , the Ricci Tensor R^ = R p wv 
and its trace R = g a/3 R a /3, namely the Ricci scalar. 

We look for static, (pseudo-) spherically symmetric (SSS) solutions with various topologies, and 
write the metric element as 

dr 2 o f dp 2 2 



ds z 



-e 2a{r) B(r)dt 2 



(3) 



B{r) \l~kp 2 

where a(r) and B(r) are functions of the radius r, and the manifold will be either a sphere ^2, 
a torus Ti or a compact hyperbolic manifold Y2, according to whether k = 1,0,-1, respectively. 
With this Ansatz, the scalar curvature and the Gauss Bonnet read 

R = -3 



d „ , , 



r dr 



—at(r)-2B(r) 
dr 

4 — — — a(r) - 
r dr 



Tr a{r) 



^B(r)-2B(r)^a(r) 



B(r) 



2k 



(4) 
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da{r)\ ( dB{r) 



dr 



dr 



(5B(r) - 3k) 



dB{r) 
dr 



d 2 B(r) 
dr 2 



(B(r) - k) 



-2(B(r) - k)B(r) 



da(r) 
dr 



d 2 a(r) 
dr 2 



(5) 



By plugging this expressions into the action ([1} , one obtains a higher derivative Lagrangian theory. 
In order to work with a first derivatives Lagrangian system, we may use the method of Lagrangian 
multipliers (for the case of FRW space-time see Refs. |2H [22 , 23J). This method permits to consider 
as independent Lagrangian coordinates the scalar curvature R, the Gauss Bonnet invariant G and 
the quantities a(r) and B(r), appearing in the topological spherically static symmetric Ansatz. 

By introducing the Lagrangian multipliers A and \i and using Eqs. (J4j)- d5j) , the action may be 
written as 



dt / dr 



(e a ( r V 2 ) ^F(R, G)-X R + 3 Qj-fl (r)\ ij-a (, ; 



+2Du-)[±a{r)) + ^B (r) + 2 B (r) (r) 



,B{r)d . . n B(r) 2k 
r dr r z r z 



Ad_ 

r dr 



B(r) 



G - 



fda(r)\UBp; ){ . dni) :U) 



\ dr 



dr 



dB(r) 
dr 



+2(B(r) - k)B(r) 



da{r) 
dr 



d 2 a(r) 
dr 2 



Making the variation with respect to R and G, one gets 

d 



A 



dR 
8 



T(R, G) 



V = —T G (R,G), 



d 2 B(r) 
dr 2 



(B(r) - k) 



(6) 



(7) 



(8) 



Thus, by substituting this values and by making an integration by part, the total Lagrangian Jz? 
of the system takes the form 



J2f(a, da/dr, B, dB/dr, R, dR/dr, G, dG/dr) = e"^ jr 2 (T - T' R R - T' G G) 



2 /-I, I A- - r^l B(r) 



dr 



dG / dBW 
GG dr V dr 



„ dR 2 



dB(r) 
dr 



2B(r) 



da(r) 
dr 



8B(r)^-)iBir)-k) 



(9) 



Now, T{R, G) has been replaced with T and we have used the following expressions: 

P = — P = — (10) 
+R- dR > QG m 

and so on. The lagrangian depends on the first derivatives of the variables at most. It is also 
easy to see that, if P G = Const, i.e. the Gauss Bonnet simply is an additive term, its contribute 
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vanishes. Making the variation with respect to a(r) and with respect to B(r), one finally gets the 
Equations of Motion: 



r\T- F' R R-F' G G)+2F' R 



-2r 2 B(r) 



»(r) ( M r ) 
\ dr 



d 2 P R 
dr 2 



4(3B(r) - k) 



dB{r) 
dr 



dB(r) 
dr 



dr 



B{r) 



dr 



dB(r) 



dr 



ArB(r) 



8B(r)(B(r)-k)^£- 



= 0. 



dT' dT' 
jWr -4(3B(r)-k) G 



dr 



dr 



d 2 V 

r 2 u ±R 

dr 2 



•4(B(r) 



k) d 2 P G 
dr 2 



(11) 



= 0. (12) 



The above equations with Eqs. (U])-© form a system of four ordinary differential equations in 
the four unknown quantities a(r), B(r), R — R{r) and G — G{r). In principle, given the model 
J-(R,G) as a function of r, our expressions permit to derive both a(r) and B(r), that is the 
explicit form of the metric. On the other hand, it is also possible to try to reconstruct the models 
by starting from the solutions. For example, fixing the form of a(r) one may reconstruct the 
model by using Eq. (IT2^) and therefore B(r) from Eq. (fTTj) . In general, the Lagrangian F(R) one 
eventually finds is not unique since one has to infer its form starting from the value it assumes 
on the solutions (see Section 4). We observe that by explicitly written R and G as functions 
of B(r) and a(r), we reduce the system to the usual high order differential equations system of 
F(R, G)-gravity. Since in this paper we would like to explore (topological) SSS solutions in Gauss 
Bonnet gravity, in what follows we will always consider the cases of non trivial contribution of G 
into the action. 



3 Solutions with constant a(r) 

In order to see how the procedure of reconstruction works, we first divide Eq. (jlip with respect to 
e a(r) r 2 anc j th en perform the derivative with respect to r, obtaining 



dr 



R 



dr R 



3B(r) 



( my) 

\ r 2 



2k QdB{r) d 2 B(r) 



2k 1 d 2 B(r) 
r 3 r dr 2 

d 2 



d?_ 

dr 3 



dr 2 ^ G 



f28B(r) dB{r) 



dr dr 2 

3 dB(r) 6B(r 



dr 2 ^ R 



dB(r) 4B(r) 
dr r 
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dr 3 



F'c B{r) 



dr 
B(r) 



12kdB(r) 16B(r) 2 
r 2 dr r 3 

= 0. 



2k\ dB(r) 
r 3 ) dr 



VakB{r) 



3B(r) k \ d 2 B(r 



dr 2 



(13) 



The utility of this equation is due to the fact that it depends on F R and Tq only, which may be 
derived from Eq. (fT2)) ones a(r) is fixed. Let us consider the important case of a(r) = (it is 
equivalent to a(r) = Const) into the metric ©, namely 



ds 2 = -B(r)dt* + ^ F) +r> + P 2 ^ ) • ( 1 4 > 



The equations ©-© lead to 

dr 2 r dr 



d 2 , x AdB(r) B(r) 2k , , 
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_ 4 fdB(r)Y 4 d 2 B(r).. . ,. 

From Ea.(fT2j) one has 

-r 2 ^< R + 4(B(r)-k)^ = 0. (17) 

A simple choice is to investigate the models with d 2 J r n/dr 2 = and d 2 F' G /dr 2 = 0, namely 

F' R = ar + b, T' G = cr . (18) 

Here, a, b and c are integration constants (the pure constant in T' G does not give contribution to 
(fT5)) and will be analyzed in the next Section). Eq. (JT3J) leads to 




(2B(r)-k) + ^{B(r)-k) 



adB(r) 
r 



dr 



4(3S(r) - fc) dB( 




= 0. 



(19) 



Unfortunately, this equation can not be solved for general values of the parameters. Otherwise, it 
may be used to find some specific solution. For example, if a = and b = 1, a solution is 

B{r) = -k + Cr, (20) 

with c = 1/(4C). From Eqs. H5])-C!l|) one nas 

J.-**^* (21) 

.22) 

If we use the fact that r = 2C/yG, we can finally solve Eq. (fTTj) with respect to ^(i?, G). As a 
result, we find that the model 

T{R,G) = R + VG , (23) 



generates the (topological) SSS solution (fT4")) with J3(r) given by ([20")) . 

4 Solutions with constant Gauss Bonnet invariant 

We will check for topological SSS solutions wich lead to a constant value of the Gauss Bonnet. We 
start by considering the case a(r) = of the metric (fl"4")l . By putting G = Go, where Go = Const, 
in Eq. one derives the form of B(r), 

B(r)=k±^ + C )±^, (24) 

where C\ and Ci are integration constants. Now, we have two very simple cases. The first one 
corresponds to T' R = 0, it means that the model J-(R, G) — F(G) depends on G only and we deal 
with a pure Gauss Bonnet theory. Thus, Eq. (|T2"j) trivially is satisfied and Eq. (ITT1) reads 

F(G o )-G o i^(Go) = 0. (25) 
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Of course, there are infinity functions F(G) which reduces (on shell) to Eq. (j2l)]) for G = Go- For 
example, one finds that the specific model 



F{G) = 7 e G/Go , (26) 

where 7 is a dimensional constant, exibits the topological SSS solution (jl4|) described by B(r) in 
Eq. ([2~i]) . but also the class of models 

F(G) = 7 G" + A , (27) 

where A = 7(71 — 1)Gq, exhibits the same kind of solution. 

An other simple possibility is to take constant the Ricci scalar also, R — Ro, where Rq — Const. 
As a consequence of Eq. dU, it follows that C\ = C% = in the solution (|M|) and Ro — y/6Go- 
Also in this case, Eq. (fT2")) trivially is satisfied and Eq. (TTTJ) reads 



T{Ro, Gq) - RqF' r {Rq,Go) — G Q J r ' G (RQGo) + 2J r ' R (R ,Go)J — . (28) 
It means that in principle every F(R, G)-model admits the de Sitter solution 

B(r) = (k-^p), (29) 



where A = ^SGo/8 is fixed by Eq. ([2"g]). such that R = 4A and G = 8A 2 /3. For example, the 
model F(R : G) = R + jG 2 , with 7 a dimensional constant, exhibits the de Sitter solution with 
A = [9/(327)] 1 / 3 . This is a well-known result. 

Finally, we observe an interesting invariance of the EOM, occurring when Go = 0. If the model 
assumes the form 

F{R, G) — Gf(R, G) , (30) 

where f(R,G) is a function of R and G again and limc^o Gf(R, G) = 0, the two Equations of 
Motion (fTTjl - (|12p trivially are satisfied. Let us have a look for some examples. 

Case a(r) =0 and Go = 

If a (r) = 0, by putting G = in Eq. (|2"3)l. we get 



B(r) = fc±Civ/l + G 2 r, (31) 

where we have redefined the constants. The models in the form of (f3"0"|) admit this kind of SSS 
solutions. 



Case a(r) ^ and G = 

We simply can also verify that Go = provided by 

B(r)=k, k = l, (32) 

independently on a(r). In the topological case k = 1 (the sphere), the models in the form of (|30p 
admit every SSS solution with B(r) — 1 and a(r) a general function. This class of solutions can 
not describe the black holes. 



fi 



5 Black hole solutions and Wald entropy 



We may define a black hole event horizon as soon as there is positive solution th of 

dB(r) 



B(r B )=0, 



dr 



> 0. 



(33) 



The second condition needs to be imposed to preserve the metric signature out of the horizon. So, 
for the metric in Eq. (|2H)) , it turns out that only in the topological case with k = 1 (sphere) we can 
describe a black hole if C > 0. For metrics in Eqs. (f2"4")) -(|3i p . we can have black hole solutions by 
choosing the appropriate sign in B(r) according with the topology. In what follows, it is always un- 
derstood that we will suppose to deal with SSS solutions which effectively describe the black holes. 

The entropy associated to the black holes solutions can be calculated via Wald's method [18J. 
Following Refs. pHl I25| . the explicit calculation of the black hole entropy Sw is provided by the 
formula. 

ffJSP 



hi dp d(j) . 



(34) 



H 



We use the suffix 'H' for all quantities evaluated on the horizon. Here, S£ = ^(R^a, R^, R, 
dfiviV-yRfivZa...) is the Lagrangian density of any general theory of gravity. The antisymmetric 
variable — —e utl is the binormal vector to the (bifurcate) horizon. It is normalized so that 
e^e'"' = —2 and, by using the metric Ansatz ([3jl. it turns out to be 



(35) 



5* being the Kronecker delta. Finally, the induced volume form on the bifurcate surface r = 
rjf , t = const is represented by \fh~2 dp d(f>. The variation of the Lagrangian density with respect to 
Rfiv^a is performed as if R^^a and the metric g^ v were independent. Since for the models under 
consideration _£f = J 7 (R,G)/(2k 2 ), formula (|34p becomes 



S w = -8nA H e 2a{rH) 

8TrA H e 2a{r ^ 
2^2 



<5-Roioi 
5R 



H 



SRo 



6G 
SR 010 



H 



(36) 



Above, Ah = Vkf H , in which V\ — Air (the sphere), Vq = |3t|, with r the Teichmiiller parameter 
for the torus, and finally V—\ — 47rg, g > 2, for the compact hyperbolic manifold with genus g [26J. 
Since 

SR 1 



SG 



= [2R^° - 2{g^R va + g™R^ - g^ R v * - g^R^) + - 9^ a g vi )R] 



SRfj,u^a 

by using the horizon condition B(ru) = one obtains 



Sw 



A 



H 



4G 



N 



H 



(37) 
(38) 

(39) 



In General Relativity T(R, G) = R and one recovers the Area Law, namely Sw = Ah /4G 



N ■ 
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6 First Law of thermodynamics 



Evaluating the first equation of motion of J-(R, G)-gravity (llip on the event horizon, and multi- 
plying both sides of the equation by Vk/(2n 2 ), we directly obtain 

e a <r») dB(r) \ dS w 



4-7T dr 



h drn 



a ( rH ) ( kT' R (RH,GH) _ RhF'r{Rh,Gh) + GhF' g (Rh, Gh) - F(R Hi Gh) 2 \ Yt (ac\\ 
6 ^ 2G N AG N Th ) An ' 1 UJ 

where Rh and Gh are the Ricci scalar and the Gauss Bonnet invariant evaluated on the horizon 
and Sw is the Wald entropy given by Eq. (|3"9")). Here, we have used the BH horizon condition 
B{r H ) = 0. 

In 1974 Hawking found an important correlation between the proprieties of black holes and the 
laws of thermodynamics [20]. The black holes can irradiate due to quantum effects, that is the so 
called Hawking radiation. For static black hole described by the metric the Killing/Hawking 
temperature of the related horizon reads 



e a M dB(r) 



An dr 



,41) 



This is a well-known result, and it can be justified in several ways, for example making use of 
standard derivations of Hawking radiation |27| , or by eliminating the conical singularity in the 
corresponding Euclidean metric, or either by making use of tunneling methods, recently introduced 
in Refs. [28 , 29], and discussed in details in several papers. 

Thus, in analogy with the case of F(i?)-gravity treated in Ref. [19] , one may try to derive the 
First Law of T{R, G)-black hole thermodynamics from Eq. (|4T))) . where the Killing temperature 
appears in a natural way. If the entropy depends on rjj only (it means, when T' R G ^ Const, Rh 
and Gh do not depend on the integration constants of the solution), we can write 

TkASw — 

a{r H ) ( ^-^rKRh-, Gh) RhJ~' r (Rh,Gh) + GhJ~' g (Rh,Gh) — F(Rh,Gh) 2 \ ^4 j / A0 \ 
{ 2G^ 4G^ rH )^ drH - {A2 > 

In this case, the First Law of J-(R, G)-black hole thermodynamics reads 

AE K := T K AS W , (43) 

where AEk is the variation of Killing energy Ek , and 

E K : = 



a a(r ff ) f k J~' r (Rh, Gh) _ RhJ~' r (Rh,Gh) + GhJ~ g (Rh,Gh) F(Rh,Gh) 
4ttJ \ 2G N AG N 

Of course, all this is reasonable only if th depends on a unique variable parameter, which will be 
identified with the mass of the black hole. If rjj and as a consequence also the entropy depends on 
other variables, then other thermodynamics potentials will appear and the expression of energy 
will be different. 

Unlike to the case of pure F(i?)-gravity [T7], in Gauss Bonnet modified gravity the BH entropy 
may depend on the integration constant of the solution. It is the case of the model (|23l) , for which 
the first law of thermodynamic can not been recovered by starting from the EOM, being Eq. (|44l) 
not valid. However, since in this case the parameter C — C{th) which appears in (f2"0")) is the 
solution of B(th) = 0, we may evaluate the variation of entropy as 

dC 

AS w (r H ,C) = d rH S w (C,r H ) + dcS w (C,r H )-^ , (45) 
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and the First Law of thermodynamic holds true as soon as Eq. ([45]) can be written. In this 
example, we will see that the First Law permits to identify the mass with a quantity proportional 
to the constant of integration, which explicitly appears in the black hole solution. 



Examples 

Let us consider the model F{R, G)=R+ \fG of Eq.pP]) with the SSS solution B{r) = -1 + Cr 
and G = 4C 2 /r 2 , which corresponds to the topological case k = 1 of ([120]) . We remember that for 
k = 0, — 1 the metric does not describe a black hole. The BH event horizon is defined by 

TH = £ , (46) 
such that C = 1/rjj- In this case the First Law leads to 

±E K :=pL, (47) 



and one has 



permitting the identification of the energy with the integration constant of the solution. 

In the case of the models with constant Gauss Bonnet solution, we have from Eq. (|59"|) that 
the contribute of the Gauss-Bonnet to the entropy is a constant. As a consequence, for the model 
F{G) = 7 e G / Go or F(G) = ~/G n +A, A = j(n-l)G^, with BH solution (H]), one obtains AS W = 
and Eq. (|4"3")) gives us AEk = 0. Also for the class of models (f5U|) with solution (f3"TT) . we can use 
Eq. (|44p to evaluate the Killing energy and it is easy to see that the result is Ek = (in this case 
Sw = 0) and the Firts Law trivially is satisfied. 



7 Maxwell SSS solutions in Gauss-Bonnet modified gravity 

In this Section, we will extend our formalism by considering topological SSS J-(R, G)-solutions 
in the presence of Maxwell field. We will also provide an exact solution for a specific model of 
Gauss-Bonnet modified gravity. The Maxwell action reads 

Iem = 7 / </ ; r x ~ /•••••/;., , V k F k > =0, (49) 

where run from to 3, V M is the covariant derivative operator associated with the metric 
g^v and i*y is the electromagnetic field strength. In the case of metric ([3]), due to the symmetry, 
it is easy to see that the only non vanishing component of the electroagnetic field is 

e Q(r) Q m e- a ^Q , x 

F oi = ^2^> ^ 01 = (50) 

Q being the electric charge of electromagnetic field. The model is now described by the action 

/=^2 / d A xJ—gF{R,G) + - A [ ,/ ; , ■, —/••'••/••,... (51) 

By using the relations ([50]) and by following the same lagrangian derivation of Section 2, we have 
that the EOM are equals to (Tn]) - ([T^]) with an additive term in (TTTj) due to the electromagnetic 
field: 
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*(r) 



r 2 {F - F' R R - T' G G) + 2T' R 



k- 



fdB(r 



-2r 2 B{r) d -^+A{ZB{r)-k) 
dr z 



(r) ( da ( r ) 
\ dr 



2rT' 



\ dr 
dB(r)\ dP G 



-B(r) 



dr 



dr 



dr +mr)(B(r)-k) ^ 



dB(r) 
dr 

d?p G g 2 



dV 
dr 



dT' 

4(3B(r)-fc)=7^ 
dr 



d 2 V 
dr 2 



4(B(r) - fc) 



2r 2 



dr 2 



4rB(r) 
= 0, (52) 
I = 0. (53) 



Here, we have defined Q 2 = 2k 2 Q 2 . 

Let us see for an explicit example and consider the pure Gauss-Bonnet gravity model ^(R, G) = 
F(G), which depends on G only, with the metric Ansatz a(r) = of (fT^|) . Eq. ([53")) leads to 

F G =cr, (54) 

c being a constant. By dividing Eq. (|52[) with respect to 2e"^ r V 2 and then performing the derivative 
with respect to r, we obtain 

' 4 fdB(r)\ 2 , 4B(r) d?B(r) 4(3£?(r) - fc) dB(r)\ Q 2 

~i — H o r^5 5 ; — "I — r = • ( 55 1 

dr / r z dr z T" 3 dr I r° 

In the topological case A: = it is possible to find the general solution of such equation, namely 



B{r) = -^=Jc 1 -C a r*-^. (56) 
VlO V cr 



The Gauss Bonnet invariant reads 



such that we can write r as a function of G, 

:«^)"V 5 . (58) 

and by using Eq. (|52|) we finally obtain that the model 

F(G) = - 1 -(^f-6cC 2 y°Q 2 /\ (59) 

exhibits the exact SSS Maxwell solution p^|) described by (|56p for the topological case fc = 0. 
Note that Ci is a fixed parameter of the model and G\ is the free constant of the solution. 



8 Conclusions 

If some extended theory of gravity lies behind our universe, the investigation of mathematical 
structure of this kind of theories becomes of great interest. In this paper, the lagrangian derivation 
of the Equations of Motion for studying topological SSS metrics in JF(R, G)-gravity has been 
presented. Due to a suitable form of the EOM, the study of the solutions results to be sensibly 
simplified. In particular, we found a new non trivial solution which may also describe the black 
holes in the topological case fc — 1. In this case, we defined the BH mass with a quantity 
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proportional to the constant of integration, which appears in the explicit solution, making use of 
the First Law of black hole thermodynamics and evaluating independently the entropy via Wald 
method and the Hawking temperature via quantum mechanical methods in curved space-times. 
We also showed that in principle may be possible to recover the First Law from the EOM, but, 
differently to the case of F(i?)-gravity, the validity of this derivation is restricted to few (trivial) 
cases. Finally, an exact solution for a specific Gauss-Bonnet modified gravity model in the presence 
of Maxwell field has been found. 
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